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Abstract. As it was argued by Anderson [Science 177, 393 (1972)], the 'reductionist' hypothesis does not by any means 
imply a 'constructionist' one. Hence, in general, the behavior of large and complex aggregates of elementary components 
can not be understood nor extrapolated from the properties of a few components. Following this insight, we have simulated 
different 'aggregates' of logistic maps according to a particular coupling scheme. All these aggregates show a similar pattern 
of dynamical properties, concretely a bistable behavior, that is also found in a network of many units of the same type, 
independently of the number of components and of the interconnection topology. A qualitative relationship with brain-like 
systems is suggested. 
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INTRODUCTION 

One of the most challenging scientific problems today is to understand how the millions of neurons of our brain give 
rise to the emergent property of thinking [ 1 ] . Different aspects of neurocomputation take contact on this problem: how 
brain stores information and how brain processes it to take decisions or to create new information. Other universal 
properties of this system are more evident. One of them is the existence of a regular daily behavior: the sleep-wake 
cycle |2|, |3|]. The internal circadian rhythm is closely synchronized with the cycle of sun light. Roughly speaking and 
depending on the particular species, the brain is awake during the day and it is slept during the night, or vice versa. 
All mammals and birds sleep. There is no a well established law relating the size of the animal with the daily time it 
spends sleeping, but, in general, large animals tend to sleep less than small animals. Hence, at first sight, the bistable 
sleep-wake behavior seems not depend on the precise architecture of the brain nor on its size. If we represent the 
brain as a complex network this property would mean that this possible bistability should not depend on the topology 
(structure) nor on the number of nodes (size) of the network. 

So, on one side, it has been recently argued in [4] that the distribution of functional connections p(k) in the human 
brain, where p(k) represents the probability of finding an element with k connections to other elements of the network, 
follows the same distribution of a scale-free network |5[] . Thus, in that work [ 4] , the human brain is divided in 34x64x64 
sites (called voxels) and the magnetic resonance activity of all the voxels is recorded. Let us observe at this point that 
if the human brain has around 10 11 neurons, then a voxel has around 10 s neurons. The calculation of the correlation 
matrix among the set of the voxels activities shows a power law behavior, p(k) ~ k~Y, with y around 2. This finding 
means that there are regions in the brain that participate in a large number of tasks while most of the regions are only 
involved in a tiny fraction of the brain's activities. 

On the other side, it has been shown by Kuhn et al. [6] the nonlinear processing of synaptic inputs in cortical 
neurons. They studied the response of a model neuron with a simultaneous increase of excitation and inhibition. They 
found that the firing rate of the model neuron first increases, reaches a maximum, and then decreases at higher input 
rates. Functionally, this means that the firing rate, commonly assumed to be the carrier of information in the brain, 
is a non-monotonic function of balanced input. These findings do not depend on details of the model and, hence, are 
relevant to cells of other cortical areas as well. 

Putting together all these facts, we arrive to the central question that we want to bring to the reader: are we able 
to reproduce the bistability in a complex network independently of the topology and of the number of nodes?. The 
answer is 'yes' (Fig.[TJ>. What kind of local dynamics and coupling among nodes must be implemented in order to get 
this behavior?. In the next section, we give a possible strategy for the coupling and the local dynamics which should 
be implemented in a few or many units network in order to find bistable behavior. An example of this implementation 
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FIGURE 1. QUESTION: Is it possible to implement some kind of coupling and nonlinear dynamics in each node of a complex 
network in order to get Instability?. ANSWER: Yes. 



is given in [7:]. Here we center our attention in the bistability present in the case of a few coupled functional units. Four 
different discrete models in two and three dimensions are collected. In view of the results, we want to suggest with all 
these examples the possibility of 'constructivism' in the world of complex systems. 



MODELS OF A FEW COUPLED FUNCTIONAL UNITS 



General model 



Our approach consider the so called functional unit, i.e. a neuron or group of neurons (voxels), as a discrete nonlinear 
oscillator with two possible states: active (meaning one type of activity) or not (meaning other type of activity). Hence, 
in this naive vision of the brain as a networked system, if x\, with < x' n < 1, represents a measurement of the ith 
functional unit activity at time n, it can be reasonable to take the most elemental local nonlinearity, for instance, a 
logistic evolution [8], which presents a quadratic term, as a first toy-model for the local neuronal activity: 
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FIGURE 2. Discrete nonlinear model for the local evolution of a functional unit. 
It presents only one stable state for each Then, there is no bistability in the basic component of our models. 
For pi < 1, the dynamics dissipates to zero, x'„ = 0, then it can represent the functional unit with no activity. For 
1 < pi < 4, the dynamics is non null and it would represent an active functional unit. This local transition is controlled 
by the parameter pj. The functional dependence of this local coupling on the neighbor states is essential in order to 
get a good brain-like behavior (i.e., as far as the bistability of the sleep-wake cycle is concerned) of the network. As a 
first approach, we can take pi as a linear function depending on the actual mean value, X l n , of the neighboring signal 
activity and expanding the interval (1,4) in the form: 

pi = pi (3X^ + 1), (excitation coupling) (2) 
or 

Pi = Pi ( — 3X^ + 4), (inhibition coupling) (3) 

with 

1 N ' 

*» = xrl>»- (4) 

Ni is the number of neighbors of the ith functional unit, and pi, which gives us an idea of the interaction of the 
functional unit with its first-neighbor functional units, is the control parameter. This parameter runs in the range 
< pi < p,nax, where p max >; 1. When pi = p for all i, the dynamical behavior of these networks with the excitation 
type coupling [7] presents an attractive global null configuration that has been identified as the turned off state of 
the network. Also they show a completely synchronized non-null stable configuration that represents the turned on 
state of the network. Moreover, a robust bistability between these two perfect synchronized states is found in that 
particular model (see lO] for more details). For different models with a few coupled functional units we sketch in the 



next subsections the regions where they present a bistable behavior. The details of the complete unfolding of these 
dynamical systems can be found in the references Ir9l IToL Hill . 



Models of two functional units 

Let us start with the simplest case of two interconnected (x n ,y n ) functional units. Three different combinations of 
couplings are possible: (excitation, excitation), (excitation, inhibit ion) and (inhibit ion, inhibit ion). 
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FIGURE 3. Two functional coupled units. 

The first two cases of coupling present bistability in different regions of the parameter space. The third case also 
shows bistability in a very narrow interval of parameter space but it requires a very fine inspection [ 12] that we shelve 
for a further work. 



Model with mutual excitation 

The dynamics of the (excitation, excitation) case |9] is given by the coupled equations: 

x n +\ = P (3y«+ l)jf„(l —x n ), (5) 
y n +i = P (3x„ + l)y„(l-y„). (6) 

The regions of the parameter space (Fig. |4]> where we can find bistability are: 

• For 0.75 < p < 0.86, the synchronized state, x+ = (x,x) = P\, with x = ^-{1 + (4 — ^) 2 }-. which arises from a 
saddle-node bifurcation for the critical value p = 0.75, is a stable turned on state. This state coexists with the 
turned off state xg = 0. The system presents now bistability and depending on the initial conditions, the final state 
can be xg or x+. Switching on the system from xg requires a level of noise in both functional units sufficient to 
render the activity on the basin of attraction of x + . On the contrary, switching off the two functional units network 
can be done, for instance, by making zero the activity of one functional unit, or by doing the coupling p lower 
than 0.75. 

• For 0.86 < p < 0.95, the active state of the network is now a period-2 oscillation, namely the period-2 cycle 
(Ps,Pe) in Fig. [4] This new dynamical state bifurcates from x+ for p = 0.86. A smaller noise is necessary to 
activate the system from xg. Making zero the activity of one functional unit continues to be a good strategy to 
turn off the network. 

• For 0.95 < p < 1, the active state acquires a new frequency and presents quasiperiodicity (the invariant closed 
curves of Fig. |4j). It is still possible to switch off the network by putting to zero one of the functional units. 
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FIGURE 4. Bistability in 2 functional units with excitation type coupling. 



Model with excitation + inhibition 

The dynamics of the (excitation, inhibit ion) case fioll is given by the coupled equations: 

x„+i = p{3>y n + \)x n (\-x n ), (7) 
y»+i = P (-3x n + 4)y n (l -y n ). (8) 
The regions of the parameter space (Pig. |5j where we can find bistability are: 

• For 1.051 < p < 1.0851, a stable period three cycle (Qi,Q2,Q3) appears in the system. It coexists with the 
fixed point P4. When p is increased, a period-doubling cascade takes place and generates successive cycles of 
higher periods 3x2". The system presents bistability. Depending on the initial conditions, both populations (x n ,y n ) 
oscillate in a periodic orbit or, alternatively, settle down in the fixed point. The borders between the two basins 
are complex. 

• For 1.0851 < p < 1.0997, an aperiodic dynamics is possible. The period-doubling cascade has finally given 
birth to an order three cyclic chaotic band(s) (A31 ^32^33). The system can now present an irregular oscillation 
besides the stable equilibrium with final fixed populations. The two basins are now fractal. 




FIGURE 5. Bistability in 2 functional units with excitation+inhibition type coupling. 



Models of three functional units 

Following the strategy given by relation (|2][3]l several models with three functional units can be established. We have 
studied in detail two of them [11] and their bistable behavior is reported here. 

Model with local mutual excitation 

Let us start with the case of three alternatively interconnected (x„,y„,z„) functional units under a mutual excitation 
scheme. 
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FIGURE 6. Three alternatively coupled functional units under the excitation scheme 

Then the dynamics of the system is given by the coupled equations: 

x n +\ = p(3y„+l)x„(l-x„), 
y n +\ = P (3Zn + l)y„(l -y„), 

Zn+l = P (3X„+1)Z„(1-Z„). 
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FIGURE 7. Bistability in 3 functional units with local excitation type coupling. 
The regions of the parameter space where we have found bistability are: 

• For 0.93310 < p < 0.95334, a big invariant closed curve (ICC) CI coexists with a period-3 orbit that bifurcates, 
first to an order 3-cyclic ICC (Fig. [7), and finally to an order-3 weakly chaotic ring (WCR) before disappearing. 

• For 0.98418 < p < 0.98763, the ICC CI coexists with another ICC C2 (see Ref. HI) that becomes chaotic, by 
following a period doubling cascade of tori, before disappearing. 

• For 1 .00360 < p < 1 .00402, the ICC CI coexists with a high period orbit that gives rise to an ICC C3. This ICC 
also becomes a chaotic band (see Ref. [ 1 1]) by following a period doubling cascade of tori before disappearing. 



Model with global mutual excitation 

We expose now the case of three globally interconnected (x n ,y n ,z n ) functional units under a mutual excitation 
scheme. 
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FIGURE 8. Three globally coupled functional units under the excitation scheme. 
Then the dynamics of the system is given by the coupled equations: 

x n+ i = p (x„+y n + z„ + l)x„(l -X„), (12) 

y n +\ = p(x„+y n +Zn + l)yn(l-yn), (13) 

Zn+1 = p(x n +y n +Zn + l)Zn(l-Z n )- (14) 




FIGURE 9. Bistability in 3 functional units with global excitation type coupling. 



For the whole range of the parameter, < p < 1 . 17, bistability is present in this system: 

• Firstly, two order-2 cyclic ICC coexist before becoming two order-2 cyclic chaotic attractors by contact bifurca- 
tions of heteroclinic type. Finally the two chaotic attractors become a single one before disappearing. 



CONCLUSIONS 

One of the more challenging problems in nonlinear science is the goal of understanding the properties of neuronal 
circuits lfl3ll . Synchrony and bistability are two important dynamical behaviors found in those circuits. In this work, 
different coupling schemes for networks with local logistic dynamics are proposed. It is observed that these types of 
couplings generate a global bistability between two different dynamical states. This property seems to be topology 
and size independent. This is a direct consequence of the local mean-field multiplicative coupling among the first- 
neighbors. If a formal and naive relationship is established between these two states and the sleep-wake states of 
a brain, respectively, one would be tempted to assert that these types of couplings in a network, regardless of its 
simplicity, give us a good qualitative model for explaining that specific bistability. Following this insight, different 
low-dimensional systems with logistic components coupled under these schemes have been presented. The regions 
where the dynamics shows bistability have been identified. Other low and high dimensional models merit a similar 
detailed inspection in the future. This study could put in evidence the possibility of 'constructivism' in the world of 
complex systems. 
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